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摘 要 : 冰川 冰 储 量 不 仅 是 冰川 的 重要 属性 ,而 且 是 核算 冰川 水 资源 及 预测 冰川 变化 的 基础 数据 ， 
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冰 储 量 及 其 变化 具有 重要 的 理论 与 现实 意义 。 目 前 冰川 储量 估算 的 主要 方法 有 


经 验 公 式 法 、 冰 厚 模型 估算 法 、 探 地 雷达 法 ;冰川 储量 相对 变化 计算 方法 有 实地 测量 法 和 遥感 监测 
法 。 通 过 系统 分 析 和 讨论 各 计算 方法 的 原理 \ 现 状 及 存在 的 问题 ,以 期 为 冰川 储量 估算 提供 方法 参 
考 。 研 究 表明 ,对 于 冰川 冰 储 量 计算 而 言 ,经 验 公式 法 适用 于 区 域 性 或 全 球 性 的 冰川 储量 估算 ; 模 
型 估算 法 适用 于 个 体 或 小 范围 冰川 储量 估算 ; 探 地 雷达 法 适用 于 人 类 易 到 达 区 域 冰川 储量 的 估算 。 
对 于 冰川 冰 储 量 相 对 变化 计算 ,实地 测量 法 适用 于 对 精度 要 求 高 且 满 足 实地 测量 条 件 的 单条 或 中 
小 型 冰川 , 遇 感 监测 法 适用 于 全 球 性 冰 储 量变 化 估算 ,但 需 改 进 算 法 和 提高 数据 空间 分 辩 率 。 目 


展 提供 了 新 契机 。 
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冰川 是 冰冻 圈 的 重要 组 成 部 分 ,被 誉 为 气候 变 
化 的 记录 器 和 预警 器 … 。 据 统计 ,全 球 冰川 和 冰 帆 
面积 (不 包括 南极 和 格陵兰 冰 盖 ) 4^] 5. 54 x 10° 
km? ,体积 为 (1.70 £0.21) x10? km (0.43 +0.06m 
SLE) ^ 。 根 据 中 国 第 二 次 冰川 编目 ,我 国 冰川 共 48 
571 条 ,面积 约 5.18 x 104 km? , 仅 次 于 加 拿 大 、 俄 罗 
斯 和 美国 ” 。 冰 川 是 我 国 极其 重要 的 固体 水 资源 ， 
我 国 西部 主要 的 大 江 大 河 都 有 冰川 融 水 补给 ,尤其 
对 于 西北 干旱 区 ,冰川 融 水 对 维系 区 域 脆弱 的 生态 
平衡 及 社会 经 济 可 持续 发 展 具 有 重要 意义 ,因此 冰 
川 研究 一 直 以 来 备 受 中 国政 府 和 学 术 界 关注 '" 。 

冰 储 量 作为 冰川 重要 的 参数 ,既是 建立 冰川 水 
文 模型 动力 模型 及 灾害 模型 的 重要 输入 参数 六 ， 
也 是 评估 冰川 变化 对 河川 径流 影响 及 制定 积极 防 灾 
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前 , 随 着 无 人 机 技术 的 逐步 应 用 ,以 及 冰川 流速 等 理论 模型 的 提出 ,这 为 冰川 冰 储 量 估 算 方 法 的 发 
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减灾 措施 的 重要 指标 ,还 是 预测 冰川 及 海平 面 未 来 
变化 的 重要 前 提 :“ 。 然 而 ,目前 只 有 极 少数 冰川 拥 
有 准确 的 厚度 和 储量 数据 ,大 多 数 冰川 冰 储 量 需 要 
通过 间接 方法 进行 估算 ,而且 冰川 研究 多 注重 冰川 
面积 变化 ,对 冰川 储量 的 研究 相对 较 少 。 为 充分 认 
识 中 国 冰川 资源 分 布 及 变化 状况 ,2013 年 中 国 科技 
部 启动 了 科技 基础 性 工作 专项 “中 国 西部 主要 冰川 
作用 中 心 冰 量 变化 调查 ”项目 ,任务 之 一 就 是 完成 
中 国 主要 冰川 作用 区 的 冰 储 量 估 算 , 认 清 中 国 冰川 
储量 的 现状 与 变化 。 近 年 来 ,国内 外 学 者 提出 估算 
冰 储 量 的 不 同方 法 ,如 刘 时 银 等 中 基于 经 验 公式 法 
估算 中 国 冰川 储量 约 (4.3 ~4.7) x 10° km? ;HUSS 
等 中 基于 冰 厚 空间 分 布 模型 估算 全 球 冰 储 量 为 
(170 x21) x10? kn? ; FREY 等 基于 经 验 公式 法 、 
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高 永 鹏 等 :冰川 冰 储 


冰川 表面 坡度 估算 法 及 冰 厚 空间 分 布 模型 估算 整个 
喜马拉雅 山 一 喀 喇 昆仑 山地 区 冰川 储量 介 于 2 955 
~4 737 km 之 间 。 本 文 初步 总 结 国内 外 冰川 冰 储 
量 及 其 变化 估算 方法 ,并 对 比分 析 不 同方 法 的 适用 
性 ,以 期 为 准确 评估 冰川 资源 总 量 及 其 变化 ,研究 冰 
川 水 资源 的 现状 及 变化 规律 ,量化 冰川 对 海平 面 变 
化 的 贡献 等 研究 提供 方法 参考 。 


1 冰川 冰 储 量 估算 方法 


1.1 冰 储 量 一 面积 经 验 公式 法 

冰川 冰 储 量 一 面积 经 验 公 式 法 是 对 传统 野外 监 
测 方法 的 补充 ,其 原理 是 理想 状态 下 气候 冰川 流动 
和 冰川 形态 保持 持续 稳定 ,通过 对 相关 参数 实地 测 
量 的 冰川 样本 进行 回归 分 析 , 构 建 特征 参数 与 冰川 
储量 相关 关系 ,从 而 估算 无 实地 测量 数据 冰川 的 储 
HES) 。 由 于 此 方法 相对 简单 ,不 需要 气象 .水文 高 
程 等 较 难 获取 的 参数 ,而 且 随 着 冰川 编目 数据 逐步 
完善 ,冰川 面积 数据 容易 获取 , 冰 储 量 一 面积 经 验 公 
式 法 成 为 估算 冰川 储量 最 为 简便 的 方法 , 表 1 汇总 
了 目前 国内 外 经 验 公 式 。 冰 川 冰 储 量 一 面积 经 验 公 


量 计算 方法 及 发 展 趋势 
式 基 本 形式 为 : 
V=cA” (1) 
式 中 :Y 表示 冰 储 量 (km ) ;4 表示 冰川 面积 (km ) ; 
c All y 为 比例 参数 。 


式 (1) 也 可 以 转换 为 冰 厚 一 面积 比例 关系 , 公 

式 如 下 : 
H - cAP (2) 

式 中 :万 表示 平均 冰 厚 (m) ;4 表示 冰川 面积 (km ) ; 
c 和 8B 为 比例 参数 ,6=y -1. 3€ 2 汇总 了 部 分 估算 
冰川 平均 厚度 的 经 验 公式 。 
1.2 冰 厚 模型 估算 法 

为 解决 无 实测 数据 冰川 冰 储 量 估算 问题 ， 
HUSS!” FARINOTTI /'! GANTAYAT 51 等 学 者 提 
出 了 模型 估算 法 , 表 3 汇总 了 目前 用 于 估算 冰 厚 的 
主要 模型 。 其 原理 是 采用 物理 模型 模拟 冰川 物 
质 或 能 量 传输 过 程 来 估算 冰川 储量 1。 基于 假设 
条 件 的 不 同 , 冰 厚 模型 估算 法 可 分 为 以 下 4 种 。 
1.2.1 浅 冰 近似 法 IEF PATERSON 提出 的 浅 冰 
层 近似 (SIA) 理论 Bl | HAEBERLI 和 HOELZELP' 
提出 沿 冰 川 表面 中 流 线 平 均 坡 度 佑 算 冰 厚 的 方法 : 


R1 不 同 地 区 冰 储 量 估算 经 验 公 式 


Tab.1 Empirical formulas for the estimation of 


ice volume in different regions 


Z^ IS AN —P 

冰川 "i BEN pue 
V= 0.027 - A? 亚洲 中 部 地 区 ERASOV P?! 
V= 0.036 1 - Al 4% 祁连山 SHI et al. (1°! 
V= 0.0597. Al? Svalbard 地 区 MACHERET and ZHURAVLEV"!!! 
V= 0.04 - Al 冰 帽 和 冰原 CHIZHOV and KOTLYAKOV!"?! 
V= 0.037 1 < A135 Svalbard 地 区 MACHERET et al. |?! 
V= 0.03 - A136 V= k. Al 欧洲 ZHURAVLEV ^! 
V= 0.021 8 - AL!” 美国 西北 部 DRIEDGER and KENNARD/?! 
V= 0.048 + Al I6 欧洲 ZHURA VLEV 6! 
V= 0.0298 - A+”? 阿尔 泰山 KIL MACHERET et al. [1 
V= 0.0285 - A135 阿尔 卑 斯 山 CHEN and OHMURA!! 
V= k. A+% V= k. Al” 全 球 MEIER and BAHR!'?! 
V= k » AL375 V= k AL” 型 冰川 区 BAHR I 
V= 0.0276 - A+% BIKA LI BAHR et al. 2?! 
V= 0.0213 . Al 全 球 VAN de Wal and WILD?! 
V= 0.04 - AL 祁连山 和 天 山 LIU et al. [231 
V= ke 4456-290] 全 球 RADIC et al. ?*! 
V= 0.036 5 - Al V= 0.053 8 - A1” 全 球 RADIC and Hock ??! 
V= [0.024 -0.042] - AU-26- 1:555] 全 球 HUSS and FARINOTTI'! 
V= k. Alb 38-146] 全 球 ADHIKARI and MARSHALL!*6! 


WV 表示 冰 储 量 (km ) ;4 表示 面积 (km? ) ;表示 可 变 系数 


表 2 不同 区 域 冰 厚 估算 经 验 公式 
Tab.2 Empirical formulas for the estimation of 


ice thickness in different regions 


AX 适用 区 域 提出 者 
H= -11.32 453. 2149? 天 山 苏 珍 [27] 
H= 5.2 +15.44%5 Ba BLL CHENU! 
H= 34,449.55 祁连山 ERK] 
注 : 及 表示 冰川 平均 冰 厚 (m) ;4 表示 冰川 面积 (km? ) 
T 
h= ee (3) 


7 fogsin( a) 

式 中 :hy 表示 冰 厚 (m) ; r 表示 沿 中 流 线 的 平均 基底 
HUM (kPa) sf 表示 形状 要 素 ;p 为 冰川 密度 (kg 

m 7) ;g 238 7] JE RE (m. + s 7) 5a 为 沿 中 流 线 
的 表面 平均 坡度 。 其 中 ,7 的 计算 参考 Haeberli 和 
Hoelzle 提出 的 冰川 中 流 线 平均 基底 剪 切 应 力 经 验 
计算 公式 : 

0.5 +159.8AH -43.5( AH)! AH<1.6 km 

í go "n 


式 中 : 7 为 平均 基底 剪 切 应 力 (kKPa) ; AH 为 冰川 最 
高 海拔 与 最 低 海拔 的 高 程 差 (m) 。 

1.2.2 基于 质量 守恒 的 冰 厚 估算 模型 ”此 类 模型 
是 基于 质量 守恒 定 律 来 估算 冰川 厚度 ,通常 当 冰 被 
视 为 一 个 不 可 压缩 的 介质 ,并 用 相对 应 的 连续 介质 


运动 方程 来 表示 。 冰 通 量 散 度 ( Vg) 的 变化 与 冰 厚 
变化 速率 (9h/691) 和 表 观 物质 平衡 (上 ) 有 关 , 可 表示 
为 : 


Vo= Sb (5) 


通过 估算 冰 厚 变化 速率 ( ah/at) 和 表 观 物质 平 
ffi (b) 的 空间 分 布 ,从 而 量化 冰川 的 物质 转换 并 依 
据 Glen 流动 定理 2 将 冰 通 量 转化 为 冰 厚 ,其 公 


n n+2 

d o (6) 
式 中 :h 表示 每 个 海拔 带 的 冰 厚 (m) ;a 表示 冰川 表 
面 平 均 坡度 (°) ;9 表示 通过 冰川 横 截 面 宽度 的 标准 
冰 通 量 (m + a7) ;4 为 流动 速率 因子 ;n 为 流动 规 
律 指数 ;p 为 冰 密 度 (kg. m 7) sg 表示 重力 加 速度 
(m*s*), 
1.2.3 最 小 值 法 最 小 值 法 的 本 质 是 将 冰川 的 厚 
度 问 题 转化 为 最 小 值 问题 ,通过 定义 一 个 成 本 函数 
来 判定 模型 与 实测 值 之 间 的 区 别 。 通 常 ,从 实测 数 
据 中 可 获得 冰川 表面 的 数字 高 程 模型 (DEM ) 数据 ， 
并 给 冰川 基 宕 地 形 一 个 初始 的 猜测 值 ,然后 运用 
Forward 模型 预测 可 观测 的 冰川 ,最 终 依 据 模型 预测 


RI 冰 厚 估算 模型 


Tab.3 Models for the estimation of ice thickness 


模型 类 型 所 需 数据 — 
OL 4 DEM SMB Vel. ah/at 
浅 冰 近似 法 Machguth( GlabTop2 ) [1 Y N N N 冰川 vett 
Linsbauer( GlabTop ) 2? Y N N N 冰川 ok te 
RAAJglabtop2 77! Y N N N iK 
质量 守恒 法 Farinotti ( ITEM ) [36] Y N N N 冰川 KHE 
GCbedstress!37! Y (Y) N (Y) 冰川 \ 冰 帆 
Maussion( OGGM ) 95) Y N N N 冰川 ` 冰 帆 
Huss( HF - model) ?! Y N N N 冰川 ok We 
Morlighem|*?! Y Y (Y) N vk 
最 小 值 法 VanPeltLeclercq: Y Y (Y) N vk 
Brinkerhoff ~ v2 4! ! Y Y (Y) N vk 
FuerstL421 Y Y Y Y 冰川 ` 冰 帽 
冰川 流速 法 Rabatel 7?! Y Y Y N 冰 
Gantayat Y N X N YI] kg 
Gantayat — v2 78] Y N Y N 冰川 vig 
RA AJgantayat [27] Y N Y N vk 
其 他 Brinkerhoff'“! Y Y Y Y 冰 
GCneuralnetL45] Y N N N 冰 


N 表示 不 需要 该 数据 


E OL 表示 冰川 轮廓 ,DEM 表示 数字 高 程 模 型 ,SMB 表示 表面 物质 平衡 ,Vel. 表示 表面 冰 流 速率 ,0h/9t 表示 冰 厚 变化 速率 ,Y 表示 需要 该 数据 ， 
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值 与 实测 值 之 间 的 不 同 通过 重复 迭代 (满足 最 小 成 
本 了 羡 数 ) 对 模型 进行 校正 ,从 而 获得 冰川 厚度 数据 。 
Forward 模型 1 表达 式 如 下 : 


aH à a(b +H) 
üt ax L ax | *tB,x, «x < ,ti Xt «t, 
H(x,t) =0, X =X, X, f; LI, (7) 


H(x,t;) 2s,(x) -b(x), 
RP: 五 表示 冰 厚 (m) tt, APH A SE DDR 
止 时 间 (a) 50 RI s 分 别 表示 冰 床 和 冰川 表面 地 形 
(m) ,B 表示 冰川 物质 平衡 (m w. e. ) 5D 表示 扩散 
系数 。 

1.2.4 LFAMAROKRKM JF CUFFEY 
和 PATERSON "5 提出 的 层 流 流动 方程 , GAN- 
TAYAT $ 引入 冰川 流速 来 估算 冰川 厚度 的 模 
型 ,其 计算 公式 为 : 


X, «x €x, 


(8) 


WP: h 表示 冰 厚 (m);u, Au, 分 别 为 冰川 表面 和 
基底 的 速率 (m.， a) ;4 为 流动 速率 因子 (Pa ， 
s^  ) ;n 为 流动 规律 指数 ;7 为 平均 基底 剪 切 应 力 
(kPa) 。 
1.3 ”基于 探 地 雷达 测 厚 估算 冰 储 量 
探 地 雷达 (GPR ) 是 通过 发 射 和 接收 高 频率 、 短 
脉冲 电磁 波 ,根据 接收 到 的 反射 电磁 波 的 双 程 走时 、 
振幅 波形 和 频率 等 特征 来 分 析 和 推 朵 地 下 介质 的 
结构 与 性 质 特征 等 随 深度 的 变化 情况 。 由 于 冰川 与 
岩石 介 电 性 质 的 巨大 差异 ,在 探 地 雷达 图 像 资料 中 
很 容易 识别 冰 一 岩 界 面 的 位 置 ,从 而 获得 测 点 位 置 
的 冰川 厚度 。 基 于 此 原理 ,在 运用 探 地 雷达 获取 
冰 厚 数据 之 后 ,绘制 冰川 厚度 等 值 线 图 ,通过 空间 插 
值 方法 获得 冰川 厚度 空间 格 网 数据 * ,从 而 计算 冰 
川 和 储量 和 平均 厚度 。 其 公式 为 |: 
h=At(C,./2) (9) 
SUH th AR UIE (m) ;At 为 回 波 延 时 (us) ;Ci 为 冰 
内 波 速 (m + ps!) o 
探 地 雷达 具有 全 天 时 、 全 天候 和 不 受 大 气 传播 
及 气候 影响 的 特点 ,可 很 好 的 弥补 遥感 光学 影像 的 
不 足 , 且 得 到 的 冰 厚 数据 也 最 为 精确 。 但 雷达 数据 
获取 渠道 有 限 ,数据 处 理 复杂 ,而 且 由 于 冰川 区 自然 
条 件 特殊 , 冰 厚 测量 一 般 集中 于 冰川 消融 区 ,积累 区 
缺乏 实测 数据 ,因此 ,通过 空间 插值 得 到 整 条 冰川 的 


冰 厚 分 布 误 差 较 大 。 对 于 山地 冰川 , 探 地 雷达 测 厚 
只 适用 于 人 类 可 以 到 达 的 冰川 区 域 ,而 不 适 于 大 范 
围 冰川 测量 。 


2 冰川 冰 储 量 相对 变化 估算 方法 


2.1 直接 测量 法 

直接 测量 法 是 指 利用 花 杆 和 雪 坑 直接 实地 测量 
来 获得 冰川 表面 物质 收 支 状况 的 方法 ' ”i 。 该 方法 
通过 实地 测量 ,获取 冰川 相对 于 前 一 年 消融 季 末 的 
表面 高 程 变 化 ,结合 冰川 近 表 层 的 物质 密度 ,对 冰川 
物质 平衡 水 平 进行 估算 。 利 用 花 杆 和 雪 坑 测量 冰川 
物质 平衡 ,在 不 考虑 冰川 面积 变化 的 情况 下 ,冰川 物 
质 平衡 可 作为 冰川 储量 变化 的 近似 参考 值 。 
2.2 冰川 地 形 测 量 法 

冰川 地 形 测量 方法 主要 针对 单条 冰川 ,对 于 大 
范围 冰川 主要 依靠 遥感 手段 。 通 过 对 冰川 表面 高 程 
重复 测量 ,可 获得 冰川 储量 变化 信息 。 该 方法 假设 
冰 下 基 岩 形态 和 冰川 近 表 层 冰 雪 密 度 分 布 不 随时 间 
变化 ,利用 其 他 年 份 的 测 厚 资料 和 表面 高 程 数据 , 反 
演 冰 下 地 形 ,再 与 其 他 年 份 冰 面 高 程 数据 相 结合 , 计 
算 冰 川 储量 变化 。 目 前 ,国内 外 已 有 的 观测 手段 主 
要 包括 立体 摄影 测量 .干涉 合成 孔径 雷达 (InSAR ) 
测量 和 差分 GPS WE RABE! 。 虽 然 此 方法 
最 为 常用 ,但 也 存在 一 定 缺 陷 , 即 一 般 获 取 时 间 间 隔 
较 长 (从 几 年 到 几 十 年 不 等 ) 的 物质 平衡 数据 ,难以 
反映 冰川 冰 储 量 的 年 内 变化 。 
2.3 基于 多 源 地 形 数据 和 多 时 相 遥 感 
川 储量 

近年 来 ,雷达 技术 的 发 展 .激光 高 度 计 的 应 用 、 
摄影 测量 技术 与 遥感 科学 的 结合 ,使 冰川 变化 的 三 
维 研究 取得 突破 性 进展 ,因此 利用 多 源 地 形 数据 和 
多 时 相 和 遥感 影像 数据 结合 获取 冰川 变化 的 三 维 参 
数 , 成 为 当前 冰川 遥感 技术 研究 的 热点 之 一 。 运 用 
遥感 手段 估算 冰川 储量 主要 是 通过 多 源 地 形 数据 和 
多 时 相 遥 感 影像 来 提取 冰川 区 域 的 立体 信息 。 如 通 
过 生成 DEM 数据 得 到 每 条 冰川 表面 高 程 分 布 函 数 ， 
按 面 积 加 权 得 出 冰 储 量 “” ;将 两 期 DEM 数据 进行 
偏差 纠正 后 相 减 ,在 假设 冰川 面积 不 变 前 提 下 ,与 厚 
度 变化 相 乘 得 到 冰 储 量变 化 ”| ; 计算 不 同 高 程 下 
DEM 数据 的 差 值 , 按 高 程 与 厚度 变化 的 关系 计算 得 
到 总 体 厚度 变化 值 ,与 早期 面积 相 乘 得 到 冰 储 量变 
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2.4 重力 测量 法 

重力 法 测 厚 是 用 重力 仪 在 测 厚 点 处 测 出 重力 异 
常 值 ,再 推算 出 冰 体 厚度 的 方法 。20 世纪 60 年 代 
初 我 国 在 珠 峰 绒布 冰川 .巴基斯坦 巴 托 拉 冰 川 等 考 
察 中 普遍 采用 重力 法 进行 冰川 测 厚 ”。 重 力 测 量 
法 基于 测 点 相对 于 周围 非 冰 川 区 的 重力 异常 ,因此 ， 
当 冰 川 范 围 小 、 厚 度 薄 时 ,仪器 误差 将 高 于 重力 异常 
引起 数值 波动 ,造成 较 大 的 测量 误差 。 随 着 重力 卫 
Æ (GRACE) 的 成 功 发 射 , 重 力 测量 法 得 到 广泛 应 
JH. GRACE 通过 高 精度 的 K 波段 微波 测 距 系统 探 
测 两 颗 卫 星之 间 持 续 的 空间 距离 差异 ,对 地 球 重力 
场 微小 变化 进行 测算 ,获取 探测 对 象 物质 变化 。 地 
球 重力 场 变 化 主要 由 极地 冰 盖 和 高 山 冰川 冰雪 质 
量 、 陆 地 水 储量 变化 引起 ,通过 观测 重力 场 随 时 间 的 
变化 ,可 以 反 演 上 述 变 化 状况 =, GRACE 卫星 在 
两 极地 区 的 应 用 主要 是 监测 极地 冰 盖 厚度 变化 及 其 
对 海平 面 上 升 的 贡献 量 %1 ;在 高 山 冰川 区 的 应 用 则 
侧重 于 估算 冰川 消融 速率 、 转 化 为 等 效 冰 厚度 获取 
UKE EEE ee EO", (ELS, GRACE 卫星 只 能 提 
供 优 于 月 尺度 时 间 分 辨 率 的 物质 变化 数据 ,空间 分 


辨 率 较为 粗糙 ,难以 应 用 到 小 区 域 。 
3 讨论 


3.1 冰 储 量 一 面积 经 验 公式 法 的 适用 性 

选取 表 (1) 中 3 个 佑 算 全 球 冰川 储量 的 经 验 公 
式 和 4 个 估算 阿 尔 卑 斯 山 冰 川 储量 的 经 验 公式 分 别 
做 对 比分 析 。 在 选取 的 3 个 估算 全 球 冰 川 储量 公式 
中 ,Radi C 估算 结果 最 大 , Macheret 估算 结果 居中 ， 
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Van de Wal 估算 结果 最 小 ; 当 冰 川 面 积 小 于 4 km 
时 ,Radi C 估算 结果 和 Macheret 的 估算 结果 相同 
(图 la) 。 在 选取 的 4 个 估算 阿尔 插 斯 山 冰 川 储量 
经 验 公 式 中 ,Shi 估算 结果 最 大 ; CHEN and OHMU- 
RA 和 BAHR 估算 结果 居中 , 且 前 者 比 后 者 估算 结 
果 略 大 ;ZHURAVLEYV 估算 结果 最 小 ; 当 冰 川 面 积 ; 
21 km^ 时 ,ZHURAVLEV 和 CHEN and OHMURAa 
估算 结果 相等 ,为 1.78 km ;面积 小 于 21 km? 时 ,前 
者 的 估算 结果 高 于 后 者 ; 当 冰 川 面积 等 于 23. 9 
km? ,ZHURAVLEV 和 BAHR 估算 结果 相等 (图 lb) 。 

综 上 所 述 , 同 一 区 域 冰川 冰 储 量 一 面积 经 验 公 
式 估算 结果 虽 较 为 接近 ,但 也 存在 差别 ,表明 冰川 地 
形 特 征 ,物质 平衡 变化 等 要 素 的 复杂 性 。 经 验 公 式 
本 质 上 是 基于 统计 学 的 估算 方法 ,冰川 类 型 .规模 、 
地 形 .局 地 气候 等 细微 差异 都 会 对 参数 确定 产生 影 
响 , 因 此 不 同 经 验 公 式 参数 难以 达到 绝对 统一 。 此 
外 ,经 验 公式 具有 区 域 性 ,不 适用 于 特定 冰川 储量 估 
算 , 如 ARENDT 4&^*! iz fj BAHR F 144 条 冰 
川 实测 数据 提出 的 经 验 公式 估算 美国 楚 加 奇 山西 部 
20 条 冰川 的 冰 储 量 ,误差 竟 高 达 64% 。 由 于 经 验 公 
式 法 计算 简单 高 效 , 转 换 系数 及 面积 因子 的 变化 范 
围 宕 , 且 同 类 冰川 间 的 差异 较 小 ,因此 广泛 用 于 冰川 
编目 工作 以 及 区 域 冰川 变化 研究 。 
3.2 ” 冰 厚 估算 模型 的 适用 性 

浅 冰 近 似 法 基于 有 限 边界 的 冰 体 静 力 分 析 , 其 
基本 假设 是 冰川 整体 处 于 静 力 平衡 状态 ,因此 ,该 方 
法 适用 于 物质 平衡 和 运动 速度 变化 相对 较 小 的 冰 盖 
和 冰 幅 储量 及 厚度 估算 。 在 区 域 尺 度 上 ,限于 冰川 
类 型 气候、 地 形 等 的 显著 差异 , 昌 统 一 特征 参数 难 
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图 1 不 同 经 验 公 式 估 算 全 球 (a) 阿尔 插 斯 山 (b) 冰 川 储 量 
Fig.1 Different empirical formulas for estimatingglobal (a) and Alps (b) glacier volume 
( M-Macheret , V-Van de Wal, R-Radi C,S-Shi ,Z-Zhuravlev , C-Chen and Ohmura , B-Bahr ) 
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以 确定 以 及 缺少 有 关 物 质 平衡 的 区 域 性 数据 ,因此 
该 方法 在 区 域 研究 中 具有 较 大 的 局 限 性 。 对 于 山地 
冰川 ,可 变 参数 + 和 上 的 取 值 分 别 与 冰川 谷 壁 横向 
阻力 和 冰川 横 截 面 的 形状 有 关 ,实测 值 与 理论 值 
往往 有 很 大 差异 ,如 王 宁 练 等 ”| 基于 浅 冰 近似 方法 
估算 得 到 祁连山 八 一 冰川 中 部 地 区 冰 厚 为 120 m, 
实测 的 该 地 区 平均 冰 厚 为 90 m ,两 者 误差 达 30% 。 

质量 守恒 法 的 理论 基础 是 质量 守恒 定律 ,冰川 
物质 平衡 分 布 由 冰 通 量 散 度 和 冰川 表面 高 程 变化 共 
同 决定 ,然后 根据 冰 通 量 推导 出 冰 厚 。 该 方法 引入 
表 观 物质 平衡 概念 ,与 经 验 公式 法 相 比 ,其 最 大 优势 
在 于 可 以 估算 非 稳 定性 冰川 的 冰 储 量 。 表 观 物质 平 
衡 作为 冰川 实际 物质 平衡 与 冰川 厚度 变化 率 的 差 
值 ,只 需 满足 在 整 条 冰川 范围 内 积分 值 为 零 ,无 需 考 
虑 冰川 上 任意 一 处 的 厚度 变化 率 是 否 为 零 ,即使 冰 
川 处 于 非 稳定 状态 ,该 方法 依然 适用 。 但 模型 参数 
难以 确定 与 校正 限制 了 该 方法 的 广泛 应 用 ,尤其 是 
表 观 物质 平衡 梯度 与 实际 物质 平衡 梯度 的 不 一 致 ， 
以 及 不 同 气 候 条 件 .冰川 类 型 和 规模 等 都 会 造成 冰 
川 物质 平衡 的 差异 ,如 海洋 型 冰川 与 大 陆 型 冰川 存 
在 明显 的 物质 平衡 差异 ,海洋 型 冰川 的 物质 平衡 梯 
度 远大 于 大 陆 型 冰川 的 物质 平衡 梯度 ; 另 一 方面 
校正 系数 难以 确定 ,校正 系数 变化 0.1, 平 均 冰 厚 变 
化 9% 09 ,上 且 校 正 系 数 的 确定 还 与 前 切 应 力 的 近似 
分 布 . 冰 通 量 、 基 底 滑 动 速 度 、 流 动因 子 有 关 。 即 便 
如 此 ,该 方法 对 估算 个 别 或 小 区 域 冰川 储量 仍 有 很 
大 帮助 。 

最 小 值 法 着 眼 于 冰 下 地 形 及 表面 物质 平衡 变 
化 ,适用 于 特定 冰川 与 冰 盖 储量 的 估算 。 该 方法 最 
大 的 优势 是 引入 成 本 函数 ,其 包括 三 个 方面 :首先 是 
量化 模型 与 实测 表面 高 程 的 差异 .判别 基 岩 海拔 的 
空间 差异 及 设 定 冰川 边界 以 外 冰 厚 为 零 '] ,然后 通 
过 Forward 模型 将 其 转换 为 冰川 表面 与 基底 高 程 的 
冰 流 模型 。 但 是 ,该 方法 只 对 处 于 稳定 状态 下 的 冰 
川 与 冰 盖 适用 ,对 于 准 静 止 态 及 瞬时 态 冰 川 与 冰 盖 
储量 的 估算 还 有 待 进一步 改进 与 验证 ;此 外 ,成 本 函 
数 在 使 用 过 程 中 ,如 厚度 值 为 负 模拟 与 观测 的 冰川 
表面 流速 不 匹配 .模拟 与 观测 的 冰川 表面 物质 平衡 
不 匹配 等 也 亟待 解决 各 。 

冰川 流速 法 是 依据 层 流 理论 ,基于 遥感 影像 计 
算 冰川 表面 和 底部 冰 流 速率 ,从 而 估算 冰川 储量 的 
方法 。 近 年 来 , 随 着 遥感 技术 的 快速 发 展 ,遥感 数据 
获取 更 加 简易 方便 ,因此 ,该 方法 适用 于 无 实测 冰川 
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冰 帽 的 冰 储 量 佑 算 。 但 由 于 流速 ,冰川 表面 坡度 IB 
状 因子 等 模型 参数 的 不 确定 性 ,加 之 冰川 不 同 区 域 
(如 粒 雪 倪 、 冰 川 末 端 ) 垂 直方 向 上 的 冰雪 密度 差异 
大 , 且 在 不 同 海拔 带 之 间 也 有 和 较 大 差异 ”| ,影响 该 
方法 的 估算 精度 及 应 用 。 如 GANTAYAT 等 .运用 
该 方法 估算 喜马拉雅 山地 区 CANGOTRI 冰川 的 储 
量 , 误 差 达 + 上 18% ;HELBING °°) Xt pe T Bo A GT 
山 Unteraargletscher 冰川 野外 观测 数据 与 数值 冰 流 
模型 的 模拟 结果 ,发 现 冰川 底部 滑动 速度 占 表面 速 
度 的 50% , 寿 对 底部 滑动 不 加 以 校正 往往 会 高 估 冰 
川 储量 , 需 引 入 校正 系数 加 以 修正 , 且 有 无 实测 数据 
参与 校正 系数 的 修正 对 估算 结果 影响 较 大 , 均 方 根 
误差 可 相差 20% 。 

GCneuralnet 是 基于 人 工 神 经 网 络 方法 的 冰川 
冰 厚 估算 模型 ,其 基本 假设 是 被 冰 覆 盖 的 冰 床 地 
形 与 邻近 无 冰 覆 盖 的 地 形 相似 ,然后 选择 感 兴趣 区 ， 
获取 模型 参数 ,如 表面 高 程 坡度、 坡 向 等 ,确定 神经 
网 络 体系 结构 ,对 样本 进行 反复 训练 获得 最 佳 参 数 ， 
不 断 优化 输出 层 冰 厚度 的 精度 。 由 于 假设 当前 无 冰 
履 盖 区 地 形 是 由 表面 平坦 无 起 伏 的 冰川 侵蚀 而 形 
成 ,所 以 人 工 神 经 网 络 法 更 适合 估算 表面 平缓 的 冰 
川 厚度 。 人 工 神经 网 络 法 对 输入 样本 质量 要 求 高 ， 
计算 量 大 ,不 适合 应 用 于 大 范围 冰川 储量 估算 。 同 
时 ,人 工 神经 网 络 法 更 多 地 依赖 神经 网 络 算法 ,忽略 
了 冰川 本 身 的 一 些 物理 特性 ,如 能 将 算法 与 冰川 物 
理 特性 进行 结合 使 之 参与 样本 训练 ,将 会 在 很 大 程 
度 上 提高 该 方法 的 适用 性 和 估算 精度 。BRINKER- 
HOFF 等 "提出 运用 贝 叶 斯 理论 推断 冰 床 海拔 ,其 
基本 假设 是 在 未 知 协 方差 的 情况 下 ,冰川 海拔 和 冰 
通 量 散 度 可 用 高 斯 随机 场 进行 模拟 ,因此 该 模型 适 
用 于 具有 观测 数据 的 特定 冰川 ,而 不 适用 于 全 球 性 
冰川 储量 的 精确 估算 。 
3.3 ”冰川 冰 储 量 相对 变化 估算 方法 的 适用 性 

直接 测量 法 是 最 为 直接 和 精确 的 物质 平衡 测量 
方法 ,其 能 反映 详尽 的 年 内 冰川 物质 平衡 空间 分 布 
差异 ,也 可 以 精确 地 反映 冰川 储量 的 相对 变化 情 
Ou) ,适用 于 特定 冰川 测量 与 冰 储 量 估算 。 然 而 ， 
数据 获取 过 程 中 的 测量 误差 (随机 误差 ) 和 方法 误 
差 (系统 误差 ) 均 会 增加 冰川 物质 平衡 估算 结果 误 
28 ,如花 杆 在 粒 雪 中 的 下 沉 、 测 量 点 空间 分 布 的 合理 
性 和 代表 性 及 空间 插值 方法 误差 等 ”。 此 外 ,该 方 
法 需要 重复 实地 测量 ,对 于 处 于 高 海拔 和 寒冷 偏远 
地 区 的 冰川 来 说 ,需要 付出 较 大 的 财力 和 人 力 成 本 ， 
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难以 在 较 多 冰川 上 开展 ,在 规模 更 大 冰川 上 开展 物 
apa ede 目前 ,全 球 范 围 内 只 有 大 约 

条 冰川 拥有 超过 20 a 持续 观测 数据 ,39 条 冰川 
pa 30 a 连续 数据 ,并 且 受 监测 的 冰川 主要 分 布 
在 欧洲 和 北美 等 地 区 ,亚洲 地 区 有 覆盖 较 少 ,因此 , 数 
据 代 表 性 差 ,并且 世 界 冰 川 监测 服务 处 (WGMS ) 发 
布 的 数据 缺乏 不 确定 性 分 析 和 精度 评估 ,影响 了 数 
Jet ng ur f EU 

多 源 DEM 数据 的 出 现 , 使 获得 冰川 厚度 变化 成 
yay ge’?! 。 然 而 ,由 于 不 同 来 源 DEM 数据 在 生成 
时 自身 存在 误差 ,在 DEM 数据 空间 坐标 系统 和 分 辩 
率 转换 及 匹配 等 处 理 过 程 中 也 会 出 现 误差 ,继而 对 
冰川 表面 高 程 变 化 估算 带 来 误差 ,最 终 影响 冰 储 量 
相对 变化 的 估算 精度 '“] 。 航 天 遥感 数据 在 获取 过 
程 中 ,遥感 平台 自身 姿态 稳定 性 误差 ,卫星 传感器 解 
析 力 和 获取 能 力 不 足 导致 的 数据 失真 ,大 气 辐 射 和 
干涉 对 获取 信号 的 干扰 等 , 均 会 影响 探测 对 象 在 
遥感 影像 上 的 反映 。 采 用 立体 遥感 数据 提取 DEM 
数据 过 程 中 ,软件 平台 、 地 面 控制 点 和 提取 数据 点 的 
数量 与 分 布 .空间 内 插 方 式 、DEM 数据 空间 坐标 系 
统 和 空间 分 辩 率 等 因素 的 不 同 也 会 导致 提取 DEM 
数据 产生 差异 ,这 些 非 系统 性 误差 最 终 导致 相同 研 
究 区 域内 多 源 遥 感 DEM 间 产 生 较 大 的 空间 配 准 误 
" 反映 为 DEM 数据 间 相 对 平面 扭曲 和 垂直 差异 ， 

而 影响 冰川 冰 量 估算 结果 精度 '”。 再 者 ,如 何在 
ee DEM 数据 间 求 出 冰川 厚度 变化 ,如 何 评价 其 
精度 ,也 成 为 研究 中 需要 解决 的 难题 。 虽 然 如 此 , 遥 
感 技术 仍 为 冰川 储量 的 估算 在 一 定 程度 上 提供 了 强 
有 力 的 技术 支持 。 


4 ”结论 与 展望 


冰川 储量 是 评估 全 球 海平 面 上 升 、 冰 川 变化 对 
气候 变化 啊 应 ,水 量 平衡 等 的 重要 指标 。 本 文系 统 
总 结 并 对 比分 析 了 目前 国内 外 估算 冰川 储量 及 相对 
冰 量 变化 的 方法 ,人 研究 表明 冰川 储量 估算 方法 各 有 
优 缺 点 ,经 验 公式 法 简单 快速 ,适合 估算 区 域 性 或 全 
球 性 大 尺度 上 的 冰 储 量 ,模型 估算 法 适合 估算 其 他 
TM e 以 o Mid 工作 - add 


m SE 0 直接 
测量 法 适合 估算 对 精度 要 求 高 且 满 足 实地 测量 条 件 
的 中 小 型 冰川 ;冰川 地 形 测 量 法 适合 单个 特定 冰川 
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冰 量 变化 估算 而 不 适用 于 区 域 尺度 ;基于 遥感 手段 
估算 方法 为 全 球 冰川 储量 变化 提供 了 有 力 的 技术 支 
TÉ;GRACE 卫星 覆盖 范围 广阔 同时 会 带 来 空间 分 辩 
率 低 的 问题 ,应 用 到 小 区 域 冰 川 , 需 在 球 谐 域 和 空间 
域 上 进行 改进 ,并 提高 空间 分 辩 率 。 

随 着 新 技术 .新 理论 的 出 现 , 多 学 科 的 交叉 融合 
为 冰川 储量 估算 精 度 的 提高 提供 了 新 机 遇 。 卫 星 遥 
感 技术 的 进步 ,特别 是 更 高 分 辩 率 的 雷达 和 光学 影 
fü Hg ^ ,提高 了 冰川 面积 获取 的 精确 性 ,但 表 
磺 、 积 雪 识 别 仍 是 冰川 提取 的 难题 ,未 来 可 尝试 采用 
粗糙 集 理论 及 ICESat 卫星 波形 提高 提取 精度 ,其 中 
ICESat 卫星 搭载 的 地 球 科学 激光 测 高 仪 系统 
(GLAS) 获 取 的 数据 ,被 认为 是 目前 全 球 范围 所 能 
获得 最 准确 的 地 面 高 程 数 据 “]。 在 冰川 厚度 物探 
技术 方面 ,探测 精度 不 断 提升 ,机 载 冰川 厚度 探测 雷 
达 的 逐步 应 用 也 有 望 为 区 域 性 冰川 测 厚 提供 可 
能 '“] 。 近 年 来 ,无 人 机 也 逐步 应 用 到 山地 冰川 的 测 
grep ,其 以 小 巧 灵 活 的 优势 在 山地 冰川 表面 
特征 参数 获取 和 冰川 储量 估算 方面 展现 出 极 大 的 应 
用 潜力 。 此 外 ,诸如 GRACE 重力 卫星 SAR 等 都 已 
成 为 探测 冰川 厚度 的 有 力 工具 。 

冰川 厚度 估算 模型 的 不 断 发 展 及 冰川 编目 数据 
的 不 断 完 善 为 冰川 储量 估算 带 来 新 契机 与 数据 文 
撑 。 在 冰 厚 模型 估算 方面 ,基于 浅 冰 近似 法 .质量 守 
恒 法 等 估算 冰川 储量 已 取得 丰富 成 果 。FARINOTTI 
等 对 目前 提出 的 所 有 模型 性 能 进行 了 实验 检验 ,发 
现 模 型 估算 法 与 实测 平均 冰 厚 偏差 约 为 10% + 
24% ,因此 ,如 何 提 高 模型 估算 法 的 精度 及 简化 模型 
参数 成 为 提升 模型 性 能 的 首要 任务 ;其 次 , 浅 冰 近 似 
法 中 如 何 优化 冰川 底部 剪 应 力 的 计算 方法 是 日 后 应 
用 研究 的 重点 。 随 着 雷达 干涉 测量 技术 的 发 展 ,大 
范围 快速 获取 冰川 表面 流速 成 为 可 能 ,基于 冰川 流 
速 模型 估算 冰川 冰 储 量 已 然 成 为 模型 估算 研究 的 一 
个 新 突破 口 。 随 着 冰川 编目 工作 的 持续 开展 和 不 断 
更 新 ,收录 的 冰川 信息 越 来 越 完整 准确 ,如 何 有 效 利 
用 丰富 的 冰川 编目 资料 及 其 他 辅助 数据 研发 新 的 或 
改进 现 有 经 验 公 式 是 今后 应 持续 关注 的 问题 。 
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Methods and future trend of ice volume calculation of glacier 
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Abstract: Ice volume is not only the crucial glacier property but also the foundational data in computing glacial 
water resources and predicting glacial change. Therefore ,it is of vitally theoretical and practical significance for ex- 
actly computing glacier volume and variation. In this paper ,we analyzed and discussed the principles , developments 
and problems of existing methods in order to provide the reference for the estimation of glacier volume. At present, 
major methods of estimating glacier volume include volume - area empirical formula method , method of ice thickness 
estimation models and GPR , and computational methods to estimate relative change of glacier volume include mainly 
the field measurement ,the estimation method based on remotely sensed data and gravimetric method. The volume- 
area empirical formula method was used to estimate regional or global glacier volume simply and fleetly ;the method 
of ice thickness estimation models were suitable to estimate single or small-scale glacier volume and do not need the 
input of measured data. The GPR was applied to estimate glacier volume where human beings can reach easily. The 
field measurement with the most accurate result fits to estimate glacier volume change specifically for single glacier 
as well as small or medium-sized glaciers that require the high accuracy and meet the conditions of field measure- 
ment. The method based on remotely sensed data is applicable to estimate glacier volume change globally,so it 
needs to improve the algorithm and enhance the spatial resolution of the data. GRACE has the characteristic of wide 
coverage but has a low spatial resolution ,so it is necessary to modify the spherical harmonic domain and spatial do- 
main and improve the spatial resolution when applied to small regional glaciers. With the emergence of new technol- 
ogy and theory ,the multi-disciplinary fusion provides new opportunities for the improvement of the estimation accu- 
racy of glacier volume. The development of the estimation models and continuous supplement of the glacier inventory 
data offer new chances and data support for estimation of glacier volume. 
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